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FOREWORD 


The  improvement  of  land  mobility  technology  and  the  determination  of  future 
requirements  for  the  mobility  of  combat  and  support  vehicles  of  the  Army  requires  an 
objective  evaluation  of  existing  and  conceptual  vehicles  under  various  terrain,  weather, 
and  climatic  conditions.  The  development  of  the  AMC  '74  Mobility  Model  serves  this 
purpose.  In  this  model  various  submodels  simulate  the  interactions  of  man-vehicle- 
terrain  systems  that  determine  the  overall  effectiveness  of  vehicles  in  combat  and  support 
situations.  The  off-road  performance  of  wheeled  vehicles  is  computed  in  the  areal 
module  of  the  AMC  '74  Mobility  Model  where  the  available  tractive  force  is  calculated 
from  empirically  established  relationships  on  the  assumption  that  the  motion  resistance 
is  Independent  of  the  driving  torque.  These  relations  are  available  only  for  either  purely 
frictional  or  purely  cohesive  soils. 

At  Grumman,  driven  and  towed  tire-soil  interaction  models  have  been  developed 
under  earlier  contracts  with  TACOM  (now  TAHADCOM).  These  models  are  based  on  the 
theoretical  concept  that  soil  behavior  is  characterized  by  its  fundamental  (Coulomb)  strength 
parameters  and  its  reactions  to  tire  loads  can  be  determined  by  plasticity  theory  methods. 
The  tire-soil  interaction  models  developed  on  this  basis  have  been  validated  by  test  data  on 
a variety  of  tire  sizes  and  soil  conditions;  the  driven  tire-soil  model  simulates  tire  per- 
formance realistically  by  accounting  for  the  interaction  between  driving  torque  and  motion 
resistance.  These  tire-soil  interaction  models  are  not  restricted  to  purely  frictional  or 
purely  cohesive  soils,  and  have,  therefore,  more  general  applicability  than  the  empirical 
relationships  in  the  present  Mobility  Model.  Thus,  a vehicle  model  utilizing  the  tire-soil 
models  is  called  for  to  extend  the  applicability  of  the  AMC-74  Mobility  Model  to  cohesive- 
frictional  soils.  The  vehicle  performance  model,  prepared  under  this  contract,  is  de- 
scribed in  Section  2 of  this  report.  Additionally,  a braked  tire  model  has  been  developed 
under  the  present  contract  to  make  the  computation  of  braking  performance  consistent  with 
the  veliicle  performance  model.  Development  of  the  tire-soil  interaction  model  for  braked 
conditions  is  discussed  in  Section  3 of  this  report. 

The  soil  characterization  by  its  Coulomb  strength  parameters  in  the  vehicle  per- 
formance model  allows  the  use  of  the  model  In  practically  all  types  of  soils  (Including  extra- 
terrestrial conditions).  Nevertheless,  a large  data  bank  exists  where  the  soil  properties 
are  characterized  by  cone  index  values.  To  use  this  Information  with  the  vehicle  per- 
formance model  it  is  necessary  to  convert  cone  index  values  to  Coulomb  strength 
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parameters.  The  differential  equations  of  plasticity  for  soils,  extended  for  the  axially 
symmetric  case  of  cone  penetration,  have  been  applied  to  this  problem.  A hew  method 
that  considers  "locked  in"  stress  states  in  the  soil  during  the  process  of  penetration  has 
been  developed  for  the  evaluation  of  the  variation  of  cone  penetration  resistance  with  depth. 
The  results  of  this  research  are  reported  in  Section  4. 


ABSTRACT 


A mathematical  model  of  4x4  o££-road  vehicles  has  been  developed  for  the 
estimation  of  vehicle  performance.  The  model  uses  the  pneumatic-tire  soil  interaction 
models  developed  earlier  for  driven  and  towed  tires  and  incorporates  these  as  submodels 
in  the  vehicle  performance  model.  Vehicle-soil  interactions,  such  as  redistribution  of 
axle  weights,  due  to  the  slope  angle  and  s^plied  torque,  and  effect  of  compaction  by  the 
lead  wheel  are  taken  into  account.  The  effect  on  various  torque  transfer  mechanisms 
between  the  axles  is  also  considered.  The  computer  program  for  the  vehicle  perform- 
ance model  has  been  prepared  as  a subroutine  with  suitable  arguments  for  use  in  the  AMC 
Mobility  Model.  The  vehicle  performance  model  can  be  used  with  any  soil,  the  strength 
of  which  can  be  characterized  by  its  Coulomb  strength  parameters. 

A braked  tire-soil  Interaction  model  has  also  been  developed  for  the  estimation 
of  the  braking  force  that  the  vehicle  can  develop  under  various  soil  conditions. 

A new  method  of  analysis  of  the  variation  of  cone  penetration  resistance  with 
depth  has  been  developed.  In  this  method  incremental  penetration  is  analyzed  by  assum- 
ing that  the  stress  state  in  the  soil  produced  by  the  previous  increment  remains  "locked 
in."  Cone  penetration  resistance  profiles  can  be  converted  to  Coulomb  strength  param- 
eters by  this  method  using  a trial  and  eiTor  procedure. 
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1.  SCOPE  OF  WORK 


The  scope  of  work  as  described  In  the  present  supplemental  contract  agreement 
Includes  the  following  items : 

o Development  of  a 4 x 4 vehicle  performance  model  using  the  driven  and  towed 
tire  models  developed  at  Grumman  xmder  earlier  contracts.  The  vehicle 
performance  model  should  be  suitable  as  an  alternate  submodel  of  the  AMC  '74 
Mobility  Model. 

o Analysis  of  braking  conditions  in  soft  soils  and  development  of  braked  tire 
model  to  be  used  as  an  alternate  of  the  total  braking  force  computation 
routine  in  the  AMC  '74  Mobility  Model. 

o Development  of  relationships  for  t &rious  types  of  soils  between  fundamental 
(Coulomb)  soil  strength  parameters  and  the  soft  soil  strength  data  required 
as  input  to  the  AMC  '74  Mobility  Model. 


2.  THE  4x4  VEHICLE  PERFORMANCE  MODEL 


CONCEPT 

For  some  time  it  has  been  recognized  that  the  performance  of  single  tires  or  rigid 
wheels  cannot  be  superimposed  to  obtain  that  of  a wheeled  vehicle,  since  there  are  inter- 
actions among  the  individual  running  gears  and  the  vehicle.  Instruments  of  these  interac- 
tions are  the  soil,  the  power  train  of  the  vehicle,  and  its  suspension  system.  This  latter 
affects  vehicle  performance  primarily  through  its  controlling  influence  on  ride  dynamics. 
Since  in  the  AMC  '74  Mobility  Model  (Ref.  1)  the  effect  of  ride  dynamics  on  vehicle  per- 
formance is  considered  in  a separate  module,  interactions  of  the  suspension  system  with 
the  running  gear  assembly  have  been  disregarded  in  the  development  of  the  4x4  vehicle 
performance  model.  However,  there  are  interactions  between  the  two  axles  of  a 4 x 4 ve- 
hicle that  are  independent  of  the  suspension.  These  have  been  taken  into  account  in  the 
model  development  in  the  following  way.  The  load  distribution  between  the  front  and  rear 
axle,  defined  by  the  location  of  the  center  of  gravity  of  the  vehicle  and  force  equilibrium 
conditions,  changes  with  the  magnitude  of  grade  and  the  applied  driving  torque.  In  the  ve- 
hicle model  this  interactive  distribution  of  axle  loads  is  taken  into  account  by  computing  the 
actual  load  on  each  axle  from  the  equilibrium  conditions  for  each  value  of  the  applied  driv- 
ing torque  and  slope  angle. 

In  four  wheel  drive  vehicles  the  type  of  torque  transfer  mechanisms  affect  the 
torque  distribution  between  the  axles.  A torque  transfer  coefficient  Incorporated  in  the 
model  allows  the  consideration  of  an  interaxle  differential  (front  and  rear  axle  torques 
are  equal),  no  interaxle  differential  (front  and  rear  wheels  turn  with  same  speed, 
respective  torques  are  different  due  to  differences  in  axle  loads),  or  any  torque  biased 
transfer  mechanism. 

The  soil  is  also  an  instrument  of  interaction  between  the  front  and  rear  axle.  The 
leading  wheel  compacts  the  soil  and  thereby  changes  the  properties  of  soil  that  control 
the  performance  of  the  trailing  wheel.  This  effect  is  discussed  in  more  detail  in 
Section  4. 

The  driven  tire-soil  interaction  model  developed  under  an  earlier  contract 
computes  tire  performance  for  a given  slip  value.  In  the  AMC,' 74  Mobility  Model  the 
"slip  modified  tractive  effort"  of  a vehicle  is  computed  from  the  tractive  force  that  the 
power  train  of  the  vehicle  can  provide  at  full  throttle  at  maximum,  median,  and  minimum 
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speed  in  each  gear.  To  make  the  vehicle  performance  model  compatible  with  the  AMC  '74 
Mobility  Model  it  was  necessary  to  make  certain  modifications  in  the  general  scheme  of 
computation  of  the  earlier  tire-soil  model.  This  modified  computation  scheme  provides 
for  the  computation  of  the  torques  supplied  by  the  power  train  and  the  determination  of 
the  drawbar  pull  that  this  torque  can  develop.  The  slip  associated  with  this  torque  is 
also  computed.  With  this  modification  the  computer  program  for  the  vehicle  perform- 
ance model,  prepared  as  a subroutine,  is  a self-contained  alternate  for  the  "slip  modi- 
fied tractive  effort"  calculation  in  the  areal  module  of  the  AMC  '74  Mobility  Model. 

Note  that  if  strict  adherence  to  the  structure  of  the  AMC  '74  Mobility  Model  had 
not  been  required,  another,  more  economical,  use  of  the  vehicle  performance  model 
could  have  been  chosen  for  the  determination  of  the  maximum  speed  that  the  vehicle  can 
develop  under  the  given  terrain  conditions.  In  the  AMC  '74  Mobility  Model  relatively 
simple  formulas  are  used  for  determination  of  the  slip  modified  tractive  effort,  therefore, 
its  computation  over  the  whole  range  of  speed  requires  little  computer  time.  The  actual 
speed  is  then  determined  from  the  slip  modified  tractive  effort  vs  speed  curve  as  the  max- 
imum at  which  the  tractive  effort  required  under  the  given  conditions  is  available.  The 
4x4  vehicle  performance  model  uses  more  elaborate  theory  and  computer  techniques  for 
determination  of  the  tractive  effort,  therefore,  savings  in  computer  time  could  material- 
ize with  a different  sequence  of  computations.  In  this  sequence,  the  external  resistances 
in  crossing  a terrain  unit  would  be  determined  first  and  the  vehicle  performance  model 
would  be  used  only  for  the  computation  of  torque  necessary  to  develop  an  equivalent  draw- 
bar pull.  The  maximum  speed  of  the  vehicle  would  be  determined  from  the  computed 
torque  and  power  train  data. 

MODEL  COMPONENTS 

The  main  components  of  the  vehicle  performance  model  are  the  pneumatic  tire-soil 
interaction  models  developed  for  driven  and  towed  tires  under  previous  TACOM  contracts. 
Many  computation  schemes  in  these  two  tire  models  are  common  to  both,  therefore,  the 
computer  programs  for  these  two  models  were  integrated  in  a single  computer  program 
for  use  in  the  vehicle  performance  model.  In  connection  with  this  integration,  it  was 
deemed  desirable  to  make  a slight  modification  in  the  computation  of  tire  perfonnances 
at  low  torque.  In  the  driven  tire-soil  interaction  model  the  interface  shear  stresses  are 
computed  on  the  assumption  that  the  interface  friction  angle,  6 , is  uniform  over  the 
interface.  This  assumption  is  a reasonable  one  for  relatively  high  torque  values.  In  the 
towed  tire-soil  Interaction  model  the  interface  friction  angle,  6 , is  assumed  to  decrease 
linearly  from  a 5 ^ value  at  the  entry  angle  to  zero  in  the  center  of  the  contact  area,  and 
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then  to  decrease  further  so  as  to  reach  - at  the  rear  angle.  At  low  torque  the  interface 
shear  stress  distribution  is  closer  to  that  at  zero  torque  (towed  condition)  than  the  distri- 
bution assumed  for  driven  tires.  To  avoid  inconsistencies  and  abrupt  transitions  from 
the  zero  torque  to  the  low  torque  condition,  the  interface  shear  stress  distribution  at  low 
torque  was  modified  so  as  to  ensure  a gradual  transition  from  the  low  torque  to  the  towed 
(zero  torque)  condition.  In  the  integrated  tire-soil  performance  model,  it  was  assumed 
that  at  a torque  value  corresponding  to  6 = 0.25  fi  the  zero  torque  type  distribution 
(variable  6 ) startes  to  superimpose  over  the  £ = constant  distribution.  A gradual  transi- 
tion in  the  interface  shear  stress  distribution  from  low  torque  to  zero  torque  is  thereby 
obtained.  Typical  interface  shear  stress  distributions  resulting  from  these  assumptions 
are  shown  in  Fig.  1. 

DESCRIPTION  OF  VEHICLE  PERFORMANCE  MODEL 

In  the  development  of  the  vehicle  performance  model  the  variable  names  and  other 
designations  of  the  AMC  '74  Mobility  Model  were  adhered  to.  Nevertheless,  because  of 
the  additional  capabilities  of  the  Grumman  vehicle  performance  model  it  ^ as  necessary  to 
introduce  certain  new  designations  and  input  data  information.  These  are  as  follows; 

Vehicle  Characteristics: 


ITRSF 

= 0 

No  interaxle  differential 

= 1 

Interaxle  differential 

= 2 

Torque  biased  transfer  case 

TRQFAC 

= 

1st  axle  torque/total  torque 

It  is  recommended  that  this  information  on  the  vehicle  power  train  characteristics 
be  Included  in  the  vehicle  data  sheet  and  vehicle  preprocessor  module  of  the  AMC  '74 
Mobility  Model. 

Soil  Characteristics: 


1ST  =9 


CORES  = 
PHI 

GAMMA 

SJ 

SK 


Designates  a cohesive-frictional  type  soil 
heretofore  not  included  in  the  AMC  '74 
Mobility  Model 
Cohesion  (psi) 

Friction  angle  (degree) 

Unit  weight  of  soil  (Ibs/cu  in. ) 

Parameter  in  sheai  stress-slip  equation 
Parameter  K in  shear  stress-slip  equation 


SFAC 


Parameter  Kg  in  the  estimation  of  compactlve 
effect  of  leading  wheel.  If  no  value  for  SAFC 
has  been  assigned,  a default  value  of  SFAC  = 
0.9  is  used. 


It  is  recommended  that  these  input  data  for  cohesive  frictional  soils  be  included 
in  the  primary  terrain  descriptor  module. 

The  computer  program  for  the  vehicle  performance  model  is  written  as  sub- 
routine "TIRE"  to  be  called  from  the  AMC  '74  Mobility  Model  as  an  alternate  to  the  "Slip 
modified  tractive  effort"  computation  in  the  area  submodule.  Simple  variable,  constant, 
and  array  designations  and  dimensions  are  identical  with  those  used  in  the  AMC  '74 
Mobility  Model  (except  for  the  new  terms  listed  before).  The  subroutine  TIRE  is  called 
with  arguments  in  the  following  order: 


INPUT  VALUES: 

ATF  (NG)  BTF  (NG),  CTF  (NG) 


WGHT  (i) 
DIAW  (1) 

SECTW  (i) 
TPSI  (i.j) 

DFLECT  (I,  j) 


SECTH  (1) 

VGV  (NG,  MD) 
VGV  (NG,  MN) 
VGV  (NG,  MX) 

1ST 


RCIC  (j) 


Constant  of  quadratic  fitted  to  vehicle  tractive 
effort  curve  in  gear  NG 
Weight  on  axle  i 

Outside  wheel  diameter  of  unloaded  tires 
on  axle  1 

Section  width  of  tires  on  axle  i 

Tire  inflation  pressure  ou  axle  i,  specified 

for  soil  type  j = 1ST 

Deflection  of  tire  on  axle  i,  at  pressure 

specified  for  j = 1ST 

Section  height  of  tire  on  axle  i 

Mid  range,  minimiun  and  maximum  speed, 

respectively,  in  gear  NG 

Soil  type: 

= 1 for  fine  grained  (cohesive)  soil 
= 2 for  coarse  grained  (frictional)  soil 
= 9 for  frictional-cohesive  soil 
Soil  Strength  (cone  index)  for 
j = 1 dry  season 
j = 2 normal  season 
i = 3 wet  season 
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THETA  (K) 

Slope  angle  In  radians  for  upslope 

« 

(K  = 1),  level  (K  - 2)  and  dovrnslope  (K  = 3) 

CGR 

Horizontal  distance  of  the  center  of 
gravity  from  rear  axle 

CGH 

Height  of  center  of  gravity  (for  loaded 
vehicle) 

TL 

Axle  distance 

NGR 

Number  of  transmission  gear  ratios 

TRQFAC 

See  new  designations  (Pg  4) 

ISEAS  (1) 

Indicator  for  dry  (i  = 1)  for  normal 
(i  = 2)  and  for  wet  (i  = 3)  season 

CORES 

See  new  designations  (Pg  4) 

pm 

See  new  designations 

GAMMA 

See  new  designations 

SJ 

See  new  designations 

SK 

See  new  designations 

iP(i) 

= 1 if  axle  i is  powered 
= 0 otherwise 

NTRAV 

- 1 for  traverse,  = 3 for  average  up, 
level  and  down  travel 

EFC 

Elevation  correction  factor  for  tractive 
effort 

ITRSF 

See  new  designations  (Pg  4) 

Notes  TPSI  (i,J)  and  DFLECT  (i,j)  need  not  be  specified  for  ISl'  = 9 since 

for  1ST  = 9 the  values  of  these  tire  characteristics  have  been  assumed 
as  the  average  of  the  values  for  1ST  = 1 and  2. 

OUTPUT  VALUES: 

VG  (NG,  IV)  = ^eed  In  gear  NG  modified  by  slip 

(mlnlmiun;  IV  = 1,  mid  range:  IV  •=  2, 
maximum:  IV  = 3) 

STRACT  (NG,  IV,  K)  = Slip  modified  tractive  effort  in  gear  NG  at 

minimum  (IV  = 1),  mid  range  (IV  = 1),  mid 
range  (IV  = 2),  and  maxlmxmi  (IV  = 3)  speed, 
upslopo  (K  = 1),  level  (K  = 2),  and  down- 
slope  (K  = 3) 
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FA  (NG,  K),  FB  (NG,  K) 
FC  (NG,  K) 


FORMX  (K) 


VFMAX  (K) 


= Constants  for  quadratic  fitted  to  slip 

modified  tractive  effort  vs  speed  curve  for 
gear  NG  and  slope  up  (K  = 1),  level  (K  = 2), 
and  down  (K  = 3) 

- Maximum  tractive  effort  available  In  soil 
for  slope  up  (K  = 1),  level  (K  = 2),  and 
down  (K  = 3) 

= Sipeed  at  which  maximum  tractive  effort 
Is  available 


In  the  vehicle  performance  model  the  slip  modified  tractive  effort  Is  computed 
from  the  Input  data  In  the  following  way  (Fig,  2). 

First,  the  program  checks  v/hether  Coulomb  soil  strc^ngth  parameters  are 
available.  If  not,  the  progpram  computes  the  Coulomb  soil  strength  parameters  from  cone 
Index  values  using  the  approximate  formulas  given  In  flection  4 for  cohesive  (1ST  » 1)  or 
frictional  soils  (1ST  = 2).  Then,  a fourfold  "Do  Loop"  Is  ontcrcd  for  the  computation  of 
the  slip  modified  tractive  effort  for  the  conditions  set  forth  In  the  AMC  *74  Mobility 
Model.  In  this  loop  the  tractive  force  available  from  the  drive  train  and  the  correspond- 
ing axle  torque  Is  computed  first,  then  the  axle  loads  are  determined  with  respect  to  the 
weight  redistribution  due  to  the  slope  angle  and  applied  torque. 

The  computation  of  axle  performance  is  started  with  the  axle  that  carries  the 
lesser  load.  The  reason  for  this  sequence  of  the  computations  is  that  in  the  case  of  four 
wheel  drive  and  Interaxle  differential  the  axle  torque  Is  limited  to  whichever  of  the  axle 
torques  Is  lower.  In  the  case  of  a powered  axle,  axle  performance  computations  are 
performed  by  following  the  computation  scheme  for  a driven  tire  developed  under  an 
earlier  contract  (Ref.  2)^  while  In  the  case  of  a free  rolling  axle  the  computation  follows 
the  scheme  for  towed  tires  (Ref.  3).  The  two  computation  schemes  were  Integrated  In  a 
single  program  for  their  expedient  use  In  the  vehicle  pei'formonce  model.  Details  of 
these  computation  schemes  are  given  in  Refs.  (2)  and  (3).  The  computation  scheme  for 
driven  tires  determines  the  drawbar  pull  and  torque  values  for  given  load  and  Interface 
friction  angle  or  slip.  In  the  vehicle  performance  model  the  determination  of  the  tractive 
performance  (or  drawbar  pull)  for  a given  Input  torque  is  needed  in  the  case  of  powered 
axles.  To  this  end  an  estimate  of  the  Interface  friction  angle,  6 , Is  made  and  adjusted 
In  an  iterative  scheme  until  the  computed  axle  torque  agrees  with  the  input  torque  within 
the  allowed  tolerance. 
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FRICTIONAL  SOIL 
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COMPUTE  COULOMB 
STRENGTH  PARAMETERS 
FROM  CONE  INDEX 

COMPUTE  RELATIVE  DENSITY 
FROM  CONE  INDEX  GRADIENT 
FOR  1ST  AXLE  AND  FROM  LIMIT 
PRESSURE  FOR  2ND  AXLE 
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PARAMETERS  FOR  1ST 
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VALUES  FOR  2ND  AXLE 

1 

1 



1ST  - ELSE 
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[ ENTER  DO  LOOPS  FOR  COMPUTATION  OF  SLIP  MODIFIED  TRACTIVE  EFFORT 

DO  102 
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- 1,  NTRAV  (TRAVERSE  OR  UP,  LEVEL  AND  DOWN  SLOPE) 

DO  104 
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- 1,NGR 

(FOR  EACH  GEAR) 

DO  106 
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- 1,3 

(MIN.,  MEDIAN  AND  MAX.  SPEED  IN  GEAR) 

DO  108 

IXA 

- 1,2 

(FOR  EACH  OF  THE  AXLES) 

COMPUTE  MAX.  AVAILABLE  TRACTIVE  EFFORT 
AND  CORRESPONDING  TORQUE  FROM  POWER  TRAIN  DATA 


COMPUTE  AXLE  LOADS  FOR  SLOPE  AND  TORQUE 
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F.STIMATE  INTERFACE  FRICTION  ANGLES 
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ONE  AXLE  POWERED 


COMPUTE  AXLE  PERFORMANCE 
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AXLE  BY  TOWED  TIRE  ROUTINE 


ESTIMATE  INTERFACE  FRICTION  ANGLE 
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DRIVEN  TIRE  ROUTINE 
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ESTIMATES  FOR  1ST 
AXLE  TORQUE  1 
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COMPUTED  TORQUE  EQUALS 
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COMPUTE  COEFFICIENTS  FOR  QUADRATIC  OF  SLIP.  MOD.  TRACT.  EFF.  CURVE 


Fig.  2.  Flow  Diagram  for  tha  Computation  of  Siip  Modifitri  I'ractiva  Effort 
by  thi  Vahidt  Parformanca  Modai. 
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In  £he  case  of  four  wheel  drive  the  effect  of  Interaxle  differential  Is  taken  into 
account  by  requiring  that  the  torque  on  the  axle  carrying  the  heavier  load  be  the  same  as 
for  the  axle  carrying  the  lifter  load.  The  tractive  performaiice  of  the  axle  under  the 
heavier  load  Is  computed  on  this  condition.  If  the  interaxle  differential  is  torque  biased 
then  the  ratio  of  the  torques  on  the  two  axles  is  constant. 

If  there  Is  no  Interaxle  differential,  then  the  two  axles  turn  with  the  same  speed. 
This  conditir  r.  translates  In  the  model  to  the  requirement  that  the  slip  of  the  two  axles 
and  the  Interface  friction  angle  be  the  same.  If  the  load  on  the  two  axles  Is  different,  the 
torque  on  the  two  axles  computed  on  this  condition  is  also  different.  ?'he  difference  In  the 
axle  torques  computed  In  the  model  corresponds  to  the  phenomenon  of  torque  windup. 

The  final  step  in  the  performance  computations  Is  the  determination  of  the  con- 
stants in  the  quadratic  fitting  the  slip  modified  tractive  effort  curve.  This  Information 
is  transmitted  through  the  subroutine  arguments  to  the  main  program  for  further  use. 

USE  OF  THE  MODEL  FOR  PREDICTION  AND  PARAMETRIC  ANALYSIS  OF 
VEHICLE  PERFORMANCE 

The  vehicle  purformunce  model  described  previously  incorporates  many  inter- 
active features  that  have  been  generally  recognized  as  having  an  effect  on  the  performance 
of  various  off-road  vehicles,  yet  have  not  been  considered  In  other  vehicle  performance  . 
models,  primarily  beoauso  the  general  oonoopt  of  these  models  was  directed  toward 
simplicity.  One  of  the  useful  attributes  of  this  vehicle  performanoe  model  Is  its  capability 
to  analyze  and  assess  the  signlfioance  of  these  iotoractions  under  n wide  variety  of  condi- 
tions and  thereby  gain  further  insight  into  the  interrelationships  that  govern  off-road 
vehicle  ^rfonnance.  In  the  following  disoussion,  results  obtained  by  the  vehicle  per- 
formance nodel  for  a small  number  of  selected  combinations  of  the  Input  variables  are 
presented.  These  presentations  ore  intended  to  show  the  value  of  the  model  as  an  analyti- 
cal tool.  The  conclusions  that  may  be  drawn  from  the  presented  results  are  valid  for  that 
particular  sot  of  Input  data  and  should  not  be  construed  os  generally  valid.  A systematic 
large  scale  analysis,  covering  a wide  range  of  input  conditions,  is  needed  to  draw  general 
eonoluslons.  Such  on  analysis  is  outside  the  scope  of  this  work. 

As  mentioned  earlier,  the  computer  program  for  the  vehicle  performanoe  model 
has  boon  prepared  so  as  to  comply  with  the  present  structure  of  the  AMC  '74  Mobility 
Model.  In  this  model  the  tractive  force  available  In  various  gears  at  full  thi‘ottle  and  the 
associated  slip  ore  determined  and  a tractive  effort  vs  speed  curve  is  obtained.  In  the 
following  presentation  the  same  relations  obtained  by  the  vehicle  performance  model  are 
shown  for  various  conditions,  to^other  with  the  theoretical  tractive  force  that  the  engine 


Is  capable  of  developing.  All  presentations  refer  to  the  M-161,  l/4  ton  Army  utility  truck 
(commonly  known  as  the  "Jeep")  that  has  the  following  vehicle  characteristics. 


Weight  on  front  axle  1740  lbs 

Wel^t  on  rear  axle  1460  lbs 

Wheel  base  86  In, 

Height  of  center  of  gravity  (C.  G. ) 13  in. 

Distance  of  C.  G.  from  rear  axle  46  in. 

Tire  diameter  30.8  in. 

Tire  width  7.16  in. 

Tire  section  height  7.40  in. 

Inflation  pressure  16  psi 

Deflection,  front  tires  1.31  in. 

Deflection,  roar  tires  1.14  in. 


Figure  3 shows  the  predicted  four  wheel  performance  of  the  M-151  1/4  ton  truck 
in  sand  (average  Cl  = 36),  upslope  (grade  = 40%),  level,  and  downslope  (grade  = -40%). 

The  tractive  force  shown  in  the  Figure  does  not  include  the  tangential  component  of  the 
vehicle  weight,  therefore,  the  net  tractive  force  would  be  less  when  going  upslope  and  more 
when  going  downslope.  It  Is  interesting  to  note  that  the  tractive  performance  Is,  for  all 
practical  purposes,  the  same  up,  level,  and  downslope  even  though  the  weight  distribution 
between  the  axles  is  different  In  each  case.  This  apprears  to  Justify  present  methods  of 
tractive  force  computations  on  slopes  that  disregard  the  effect  of  Interaction  due  to  the  re- 
distribution of  axle  weights  on  slopes.  However,  further  systematic  analysis  Is  necessary 
to  generalize  this  conclusion.  An  average  cone  Index  of  36  corresponds  to  a fairly  com- 
pact sand  where  the  wheel  loads  of  the  M-151  are  not  critical.  It  Is  possible,  that  under 
marginal  trafflcabllity  conditions  the  axle  weight  redistribution  could  make  the  difference 
between  "go"  and  "no  go"  conditions. 

Another  interesting  feature  of  llg.  3 la  the  decline  of  the  tractive  force  with  in- 
creasing torque  liu  the  first  and  second  gears.  This  decline  is  consistent  with  the  results 
of  experiments  and  off-road  driving  experience.  The  application  of  excessive  torque  tends 
to  spin  out  the  tiro,  Increase  its  sinkage,  and  reduce  the  tractive  force  that  the  tire  con 
develop.  In  the  present  method  of  tractive  effort  calculation  In  the  AMC  '74  Mobility 
Model,  this  decllnJitlon  of  the  tractive  force  is  suppressed  to  allow  the  application  of  the 
"full  thruttle  concept"  to  the  determination  of  the  maximum  speed  that  the  vehicle  can 
develop  under  the  given  conditions.  Should  the  realistic  tractive  force  variation  with 
speed  predicted  by  the  vehicle  performance  model  cause  any  problem  with  the  subsequent 
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Fig.  3.  FradkiUd  Four  Wh«tl  Driva  Traotiva  Farformanoa  of  tha  M-161 
% ton  Avmy  Truck  in  Sand. 


calculations  In  the  AMC  '74  Mobility  Model,  it  may  be  necessary  to  assume  that  the 
tractive  force  la  constant  in  that  speed  range  where  decline  occurs. 

The  difference  between  the  tractive  force  that  the  engine  Is  capable  of  developing 
and  the  one  actually  developed  by  the  vehicle  Is  the  motion  resistance.  It  is  seen  that  In 
contrast  to  the  constant  motion  resistance  concept  adopted  in  the  AMC  '74  Mobility  Model, 
the  motion  resistance  strongly  varies  over  the  whole  range  of  speed.  However,  Tig.  3 
should  not  be  Interpreted  as  if  the  motion  resistance  generally  decreased  with  speed;  the 
motion  resistance  decreases  because  the  applied  torque  decreases  with  speed  In  the  "full 
throttle"  concept.  The  low  motion  resistance  at  high  speed  Is  associated  with  the  low 
applied  torque  In  that  speed  range.  The  same  torque  applied  at  less  than  full  throttle  In 
lower  gears  would  result  in  the  some  low  motion  resistunoe. 

The  interaction  between  applied  torque  and  motion  reslstonoe  Is  one  of  the  most 
Important  processes  that  affect  mobility.  The  vehicle  performance  model  presented  In 
this  report  simulates  this  Interaction  well,  while  the  constant  motion  resistance  concept  Ig- 
norea  this  Important  interaction  completely.  For  a given  vehicle  the  dependence  of  motion 
resistance  on  the  applied  torque  could  be  systematically  investigated  by  means  of  the  vehicle 
performiuice  model  und  approximate  relationships  between  motion  reslstunee  and  applied 
torque  could  be  established.  It  Is  recommended  that  such  on  Investigation  be  Initiated  und 
the  results  tncorpoiated  In  the  present  structure  of  the  AMC  Mobility  Model.  A major 
weakness  In  the  AMC  Mobility  Model  thus  would  be  ellmlnutud  und  the  structure  und  effi- 
ciency of  the  model  preserved. 

Figure  4 shows  the  tractive  porformonoe  of  the  M-151  i ton  truck  In  sund  with  only 
the  rear  axle  engaged.  A comparison  of  this  hlgure  with  the  previous  one  clearly  shows 
the  enormous  advantage  of  the  four  wheel  drive  In  sund. 

hlgure  5 shows  the  upslopo  (gx'ude  40%),  lovul,  und  downslopo  (gi'ude  ■■  -40%) 
tractive  perfoi'muueo  of  the  M-101  ii  ton  truck  in  cloy  (Cl  - 30),  while  in  Fig.  0 the  level 
ground  tractive  performimoo  of  the  four  wheel  und  rear  wheel  drive  In  clay  Is  compared. 

In  I'lig.  7 the  upslope  (grade  ■ 40%)  imrformuneo  of  the  M-lBl  i ton  tx*uok  is  com- 
pared with  a similar  vehicle  equipped  with  Interuxie  dlfferuntiul.  In  sund,  the  Interuxle 
differential  would  help  tr  uetlun  duvolupment  in  the  spued  nuigc  where  the  truutlve  force 
deollnoB  with  the  decrease  of  aiKJod  (rofor  to  Fig.  3).  Ajudysos  ixirformod  for  other  eon- 
dltlons  Indlento  that  the  tratJtlvo  ijerformonco  with  no  Interaxle  differential  Is  generally 
equal  to  or  better  than  that  with  lui  Interuxle  differential.  Thu  magnitude  of  the  torque 
windup  that  uuuurs,  If  there  Is  no  Interuxle  dlffcirentlul,  tiun  id  ho  he  estimated  by  the 
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17 


vehicle  performance  model.  Figure  8 shows  the  torque  windup  (in  terms  of  axle  torque) 
in  sand  up,  level,  and  down  slope.  Another  analysis  indicated  that  in  clay  and  level 
ground  the  sign  of  the  torque  differential  changes  as  the  gear  is  changed  from  2nd  to  3rd. 
Thus,  in  off-road  travel  where  local  slopes  and  gears  are  continuously  changing  there  is 
no  buildup  of  the  differential  torque,  a fact  that  has  been  known  from  experience  but  for 
which  no  analytical  explanation  has  ever  been  offered. 

Figure  9 shows  the  effect  of  soil  strength  on  tractive  performance  in  sand.  The 
Figure  is  a confirmation  of  the  widely  acknowledged  notion  that  soil  strength  is  the  most 
important  factor  in  off-road  mobility.  Finally,  Fig.  10  shows  an  example  of  up 
(grade  = 40%),  level,  and  downslope  (grade  = -40%)  tractive  performance  prediction  by 
the  vehicle  performance  model  in  cohesive-frictional  soils  {4>=  25®,  c = 100  Ibs/sq  ft). 

The  tractive  performance  on  level  ground  shown  in  the  Figure  is  sometimes  lower  than 
up  or  downslope.  The  explanation  of  this  seemingly  paradoxical  prediction  is  that  the 
actual  normal  axle  load  on  a slope  is  less  than  the  axle  load  on  level  ground,  therefore, 
the  slnkage  is  less  and  the  drawbar  pull  higher  (the  tangential  component  of  weight  is  not 
included  in  the  tractive  force  shown).  It  is  of  interest  to  note  that  the  combination  of  a 
small  cohesion  and  a relatively  low  friction  angle  is  more  advantageous  for  the  develop- 
ment of  tractive  effort  than  either  a purely  frictional  (see  Fig.  9)  or  purely  cohesive  soil 
with  high  strength. 

These  examples  are  but  a few  samples  of  tractive  performance  predictions  for 
selected  combinsUoas  of  input  variables.  Other  conditions,  such  as  partial  throttle  per- 
formance, inflation  pressure  and  tire  size  variations,  etc.  could  be  readily  analyzed  by 
the  vehicle  performance  model.  It  is  recommended  that  the  4x4  model  be  extended  to 
multiaxle  configurations  and  used,  in  addition  to  mission  analyses,  fur  general  parametric 
analyses  of  the  tractive  performance  of  off-road  vehicles. 
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Fig.  9.  Effact  of  Soil  Strength  on  Tractiva  Parformanca  in  Sand, 
M-161  % ton  Truck. 
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3.  TIRE-SOIL  INTERACTION  MODEL  FOR  BRAKED  CONDITIONS 


CONCEPT 

Braking  exerts  a negative  axle  torque  that  has  to  be  balanced  by  soil  reaction 
forces.  The  drag  (negative  drawbar  pull)  that  results  from  the  tire-soil  interaction  under 
the  action  of  the  negative  torque  is  the  total  braking  force  that  the  tire  can  develop.  The 
negative  torque  is  balanced  primarily  by  shear  stresses  at  the  tire-soil  interface.  Con- 
ceptually, tire-soil  interaction  under  braked  conditions  is  identical  with  that  under  driven 
conditions  but  for  the  sign  of  the  interface  shear  stresses.  The  opposite  direction  of  the 
shear  stresses,  nevertheless,  results  in  certain  conditions  that  have  not  been  considered 
in  the  driven  tire-soil  interaction  model.  These  are  as  follows. 

a)  In  the  driven  tire-soil  interaction  model  a "No  Go"  condition  is  indicated 
whenever  the  soil  reactions  wore  insufficient  to  carry  the  load,  even  if  the  entry  and  rear 
angles  reached  their  maximum  value.  In  the  braked  tire-soil  interaction  model  this  con- 
dition is  identified  with  the  development  of  the  maximum  braking  effort  that  the  bralclng 
system  of  the  vehicle  is  capable  of  developing. 

b)  In  the  driven  tire-soil  interaction  model  provisions  are  made  to  account  for 
situations  when,  under  certain  soil  conditions,  only  one,  rearward  directed  slip  lino  field 
develops.  This  situation  cannot  occur  under  braked  conditions  but  an  analogous  one 
develops  when  there  is  only  ono  forward  directed  slip  line  field.  However,  the  maximum 
value  of  the  central  angle  that  defines  the  position  of  the  singular  point  is  much  larger  for 
the  forward  directed  single  slip  lino  field  (it  equals  the  maximum  entry  angle),  than  for 
the  rearward  directed  one.  The  problems  caused  by  this  dil'fereneo  are  discussed  later. 

c)  In  the  driven  tire-soil  interaction  model  the  conditions,  when  the  soil  is  not 
stressed  to  the  limiting  plastic  state,  are  identified  as  "HARD  SURFACE  CONDITIONS. " 
These  conditions,  of  course,  arc  not  critical  for  tractive  performance  but  may  bo  critictU 
for  braking  performance.  Iii  the  AMC  ’74  Mobility  Model  a braking  coefficient  (XBRCOF) 
obtained  on  pavement  is  assigned  to  this  condition.  This  braking  coefficient  may  not 
always  materialize  on  hard  soil.  In  the  bi'aked  tiro-soil  Inteructlon  model  a braking  co- 
efficient for  hard  suriaco  conditions  is  defined  as  the  ratio  of  normal  stress  (assumed  to 
be  equal  to  the  limit  pressure)  to  the  shear  strength  of  soil  under  that  normal  stress. 
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DESCRIPTION  OF  BRAKED  TIRE-SOIL  INTERACTION  MODEL 


The  computer  program  for  the  braked  tire-soil  Interaction  model  is  written  as 
subroutine  ’’BRAKE"  to  be  called  from  the  AMC  ’74  Mobility  Model  as  an  alternate  to  the 
total  braking  force  submodel  by  the  following  arguments. 

INPUT  VALUES: 


GCW 

GCWB 

GCWNB 

WGHl’ 

DIAW 

SECTW 

TPSI 

DFLCT 

SECTH 

1ST 

RCIC 

GRADE 

CGR 

CGH 

TL 

ISEAS 

CORES 

PHI 

GAMMA 

IP 

IB 

NTRAV 

XBRCOF 

XBR 


= Gross  combination  weight 

= Gross  combined  weight  on  braked  axles 

= Gross  combined  weight  on  non-brakod  axles 
= See  INPUT  designations  (for  Vehicle 
Performance  Model,  Pg  6) 

M 

M 

tl 

II 

n 

M 

tl 

II 

II 

II 

It 

It 

ft 

M 

II 

I! 

II 

ft 

■=  Maximum  combination  braking  coefficient 
= Maximum  braking  effort  vehicle  can  develop 


23 


OUTPUT  VALUES: 


TBF(J)  = Total  soll/slope/vehlcl©  derived  braking 

force  up  0 = 1),  level  0 “ 2)  and  down  (J  = 3) 
slope 

BFGONO  = 1 if  vehicle  braking  is  inadequate  for  down 

slope  operation 
0 otherwise 

The  computation  scheme  for  the  determination  of  braking  force  is  similar  to  that 
described  In  Hef.  2 for  the  tire-soil  interaction  model  for  driven  tires.  Therefore, 
details  of  the  computations  that  are  Identical  with  those  described  in  Ref.  2 are  omitted 
here  and  only  the  fundamental  organi^^atlon  of  the  computations  is  discussed  in  the  follow- 
ing description  of  the  interaction  model  for  braked  tires. 

The  flow  diagram  showing  the  major  steps  in  the  computations  is  shown  in 
l<lg.  11.  llrst,  the  Coulomb  strength  parameters  are  computed  from  the  cone  index 
values,  if  they  are  not  given  as  input  values.  This  computation  is  osseutlully  the  same 
as  in  the  vehicle  performance  model.  A "Do  Loop"  is  entered  for  the  uemputaliun  uf 
the  braking  force  up,  level,  and  down  slope  (if  required)  and  for  each  axle.  The  weii^it 
on  each  axle  is  computed  taking  into  account  the  redistribution  of  weight  due  to  the  slope 
and  applied  braking  torque.  On  the  basis  of  the  maximum  braking  torque  that  the  brakes 
can  supply  on  estimate  is  made  of  the  interface  friction  angle  - 6 . The  slip  line  field 
computation  routine  (essentially  the  same  as  in  the  driven  tire  medol)  is  entered  and 
iteration  Is  performed  on  the  entry  angle,  a until  the  normal  stress  at  a , (Fig.  12) 
matches  the  limit  pressure,  p^.  Then  the  roar  field  is  computed  and  angle  o'^,  at 
which  the  normal  stress  matches  the  limit  pressure,  is  determined.  If  the  normal 
stress  q^  at  a j.  for  an  infinitesimal  slip  lino  field  is  higher  than  the  limit  pressure, 
then  there  is  only  one  forward  slip  lino  field  and  the  computations  aru  repeated  for  u 
slip  line  field  extending  from  The  interface  stresses  obtained  from  the  slip 

line  field  computations  are  Integrated  and  the  load,  drawbar  pull,  imd  torque  determined. 
Iterations  on  a ^ and  a ^ are  perfot'med  until  the  computed  load  agrees  with  the  Input 
load  within  the  allowed  tolerances.  If  the  loud  computed  for  the  maximum  values  of 
and  o y is  less  than  the  input  load,  a "no  go"  situation  exists  and  the  available  braking 
force  is  equated  with  the  maximum  that  the  vehicle  cun  develop.  Othoi’wlse  the  available 
broking  force  equals  the  tangential  force  comix>nent  obtained  from  the  integration  of 
interface  stresses. 
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ACCliPT  INPUT  VALUES 


IF  NOT  AVAILABLE  AS  INPUT,  COMPUTE 
COULOMB  STRENGTH  PARAMETERS  FROM  CONE 
INDEX  VALUES 


ENTER  DO  LOOP  FOR  COMPUTATION  OF 
TOTAL  BRAKING  FORCE  AVAILABLE  ON  SLOPE 
DO  IQR AD  - 1 ,3  (1  - UP,  2 - LEVEL,  3 - DOWNSLOPE) 
DO  lAX  - 1,2  (AXLE  NUMBER  1 - FRONT  2 - HEAR) 


I 


COMPUTE  MAX.  TORQUE  BRAKING  SYSTEM  CAN  DEVELOP 


I 


COMPUTE  AXLE  LOAD 


X 


ESTIMATE  INTERFACE  ANGLE  -/  FROM  TORQUE 
^ tlWH  I 


ASSUME  U,,  U^, 


COMPUTE  FRONT  SLIPLINE  FIELD 


I 


I 


ADJUST  aJ 


YES 


COMPUTE  Mr  FOH  INFINITESIMAL 
HEAR  FIELD 


Mf  <P« 


X 


COMPUTE  HEAHFIELD 
DETERMINE  aj 


Mr  ^ P» 

zxz: 


NO  REAR 
FIELD 

□c 


EXTEND  FRONT  FIELD 
TO  a. 


^li^urL^LoioTTORisoi!: 

UHAKINGTOI^E  FOR  AXLE 


■UliA 


COMPUTED  LOAD  INPUT  LOAD 


X 


ONTINUE  UO  LOOP 


X 


r 


ADJUST  Uf,  Uj 


END  OF  LOOP 

-- z.x:-: 


COMPUTE  TOTAL  VEHICLE 
BRAKING  FORCE 


Fig.  11.  Flow  Diagram  for  tha  Computation  of  Total  Vahida  Braking  Foroa 
(For  Daiignationf  taa  Fig.  12). 
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PROBLEMS  ENCOUNTERED 


Unoxpooted  problems  arose  with  the  Iteration  on  a ^ aimed  at  detei^lnlng  the 
entry  angle  from  the  requirement  that  q^  match  the  limit  pressure.  This  Iteration 
scheme  was  designed  In  the  driven  tlro-soll  Interaction  model  on  the  assumption  that 
monotonlcally  Increases  with  the  entry  angle  a^.  It  was  found  that  under  certain  brulclng 
and  soil  conditions  this  assumption  Is  not  valid.  A computational  routine  was  added  to 
the  program  that  recognizes  this  situation  and  directs  the  program  to  go  to  the  sequences 
devised  for  the  condition  that  for  a ^ mux.  the  normal  stress  q^  ut  a ^ lu  less  than 
the  limit  pressure  i)|. 

Another  problem  was  onuountered  In  situations  where  conditions  called  for  a 
single  fox'ward  directed  slip  line  field.  Although  this  situation  Is  analogous  to  the  single 
roar  slip  lino  field  oondltlou  occurring  with  drivan  tires,  the  maximum  value  uf  the 
rjunlrul  angle  to  the  singular  point  In  the  slip  line  field  is  much  tdghur  In  bruited  than  in 
driven  conditions.  As  a cunscquuncc,  the  forward  dlruutod  single  slip  line  field  extends 
over  u longer  arc  length  and,  If  It  Is  shallow,  the  "J"  lines  intersect  the  luterfuue  ut  u 
very  acute  angle.  Under  these  cundUlous  the  duturmlnutlun  of  the  louatiuns  of  tlie  nudul 
points  ut  the  luturfuue  buoomuu  inaccurate  and  uurtuln  Iter  aliens  become  divergent. 

Three  methods  were  considered  to  remedy  this  situation: 

u)  Use  of  u finer  grid  that  would  multu  the  eompututluas  mure  ueeurutc 

b)  Local  refinement  of  the  computations 

e)  Local  reduction  of  thu  Intorfucu  friction  luigle  for  the  i>urposeu  of 
compulation 

l<’or  pructiuul  reasons  method  u)  wus  adopted.  It  was  fuund  that  a minor  rcductlen 
of  the  local  intciiiico  I'rlctluu  angle  augments  the  angle  of  the  Interseetiun  uf  the  hiLeiTtiue 
and  slip  lines  suffiulently  tu  eliminate  this  eompututlonul  prublem.  Method  u),  altliuugh 
simple,  would  huvu  ruqt;iirud  u slgnllieuntly  increased  core  space  fur  the  eomputatlutiH 
which  wus  nut  readily  uvullablo  on  the  liP-'UUUO  computer.  Method  u),  however,  may  be 
considered  when  the  progx  um  Is  used  with  the  CDC  U(300  computer  uvulluble  fur  the 
Mobility  Model.  Moiliod  b)  would  have  re<4uli  ed  elaborate  programming  not  Justified  by 
the  signifleunce  of  the  problem. 


4.  USE  OF  COULOMB  BTllENGTII  PARAMETEllS  OF  SOIL  AND 
CONE  INDEX  VALUES  IN  MOBILITY  J']VALUATION 


INTIIODUOITON 

In  the  vohlclo  povl'ornuuiuo  niudol  thu  soil  is  modolud  by  its  Coulomb  strungth 
ptvrumotoi's  that,  I'or  thu  purixiso  of  mobility  evaluation,  are  suitable  for  the  oharaotor- 
izatlon  of  praotloally  every  typo  of  soil.  The  use  of  those  fundamental  parameters  allows 
the  upplluation  of  soil  mocdiatiics  theories  and,  speelfieully,  thu  plasticity  theory  to  the 
problems  of  mobility  . 

On  the  other  Imid,  fur  thu  field  determination  of  soil  properties,  eone  penetration 
tests  are  extensively  used.  Diiuenniomi}  analyses  indieated  that  if  the  soil  is  either 
purely  I'rlotional  or  purely  ooliosivo  tlani  uonu  index  values  are  sulfleiont  for  the  ehurao'- 
terizatlon  of  soils  and  simple  relationsldps  between  dimeiujlunless  tire  perfurmunee 
parameters  and  ao-oalled  "numeries"  eun  be  outabliBlmd  exi)erimontally , The  experimental 
information  avuiluide  today  on  tire  ptnToi’manue  in  these  two  types  of  soil,  that  represont 
extremes  of  thti  soil  upeutrum,  is  hivaluablo  for  mobility  ovaluation. 

Ono  reason  bohbid  oentoring  tho  oxporimontal  rosoaroh  aromid  these  two  oxtrome 
tyi>os  of  soils  was  that  those  are  also  tho  soils  whoro  "no  go"  situations  most  froiiuontly 
ooour.  Today,  whon  agility  us  well  us  gonoral  mobility  is  required  of  oonibat  imd  support 
voliielos,  oritloul  situations  for  agility  imo  more  likely  to  ooour  in  the  general  class  of 
ITletioiud-eoheslvo  soils  than  in  the  extremes  of  tlie  soil  spoetrum.  Thus,  it  is  inercius- 
ingly  linportimt  to  eharaetojTvnt)  soils  by  their  finulamontal  Coulomb  stirongth  paramuturu 
rather  than  eone  index  values. 


llEljATlCyNBIlll’B  BETWEl'lN  CONE  INDl'lX  VAIdll'lS  AND  COULOMB  HTllENTIl 
I'AUAMETEUM 


In  the  AMC]  ’74  Mobility  Model  the  terrain  Is  eliaraetorlKod  by  a single  eone  Index 
value  that  represonts  the  average  eone  index  in  tlie  uppei’  six  biolms  of  soil.  Obviously, 
a single  value  emmet  be  uniquely  related  to  sovex'ul  liulependent  parameters  mid  relutlqn- 
ships  among  them  emi  only  bo  developed  if  one  of  the  Coulomb  parmueters  is  nearly  /.evo, 
as  in  tho  ease  of  either  purely  i'rletlonal  or  ijuroly  uohoslvo  soils,  hi  those  eases  the 
I'olluwlng  uppi, ’0X1110110  rolutlonshipu  have  been  developed  earlier  (lief.  2). 


2B 


Frictional  Soils 


A relationship  between  friction  angle  and  cone  liidcx  gradient  has  been  established 
indirectly  by  using  empirical  relationships  between  relative  density  and  friction  angle  ind 
estimating  the  frictional  angle  on  the  basis  of  relative  density. 

The  relative  density  of  frictional  soils  may  be  expressed  as  (lief.  2) 

- 71.1  Log  G + 11.33  ±6.8  (1) 

where  = relative  density  (%) 

G ==  cone  Index  gradient  (i)sl/ln« ) 

The  friction  angle  ( 0 ) is  estimated  from  the  following  relationship 

cot  0 = 1.64  - 0.68  * (2) 

Cohesive  Soils 


Thu  following  ompirioal  relationship  is  used  to  estimate  the  Coulomb  strength 
parameters  for  cohesive  soils  (lief.  2). 


C (iisi)  » Cl/12. 8 
0 (degrees)  « Cl/4 


(3) 


rhe  friction  angle  of  purely  oohoslvc  soils  is  theorotioally  vioro.  However > real 
soils  that  full  in  the  category  of  cohuslvu  aullo  oxlilbit  a small  friction  angle  as  Indicated 
by  Eq,  .3.  The  use  of  a small  friction  angle  is  also  advantageous  from  the  computational 
ix>int  of  view.  A number  of  oomputatlonul  schemes  are  not  upplioable  for  a 0 ■=  0 condition 
and  the  additional  algoritlims  needed  to  provide  for  this  contlugenoy  would  lucx'oaso  the 
length  of  the  program  appreciably . 


APPLICATION  OF  PLASTICITY  THEORY  TO  THE  CONE  PENETRATION  PROBLEM 


iTastlclty  theory  may  bo  applied  to  the  determination  of  cone  penetration  rosiotimco 
in  soils  that  oxliibit  rolutlvoly  small  vohmo  ol'aiige  necessary  to  develop  tlieir  shear 
strength.  In  such  sells  the  oono  index  (CT  « cone  penoli’ation  roBlstatxce/lxaao  Jirea)  may 
bo  expressed  as 

Cl  =f  (0,  0 , 7 , 6 ,t)  (4) 

where  e cohesion 

0 « friction  angle 
7 = unit  weight  of  soil 
& interface  friction  angle 
t ~ depth  of  base  level  beneath  the  surface 
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In  the  field  determination  of  soil  properties  the  unit  wei^t  of  soil  ( 7 ) is  rarely 
measured.  If  estimated  values  of  7 are  used  a source  of  inaccuracy  is  Introduced  in  the 
above  relationship.  Another  source  of  inaccuracy  is  the  value  of  interface  friction 
angle  6 , that  depends  on  the  friction  developed  at  the  cone  face.  Very  little  is  known 
a'xiut  the  actual  values  of  6 that  are  likely  to  vary  with  soil  conditions  and  the  rougluiess 
of  the  cone  face.  The  uncertainties  In  the  value  of  7 and  6 set  a certain  limitation  to 
the  accuracy  of  any  relationship  that  can  be  developed  on  a theoretical  basis  between  the 
Coulomb  strength  parameters  and  cone  index  values.  It  is  for  this  reason  that  no  attempt 
was  made  to  refine  the  approximate  relationships  (Eqs.  1,  2,  and  3)  established 
empirically. 

The  cone  index  value  also  depends  or?  the  depth  of  base  level  beneath  the  surface. 

If  an  average  Cl  Is  given  this  depth  Is  assvjnod  as  3 inches.  In  an  earlier  Grumman 
research  project  plasticity  theory  methods  were  applied  to  the  problem  of  cone  penetration 
(Ref.  4).  An  interactive  computer  prograin  was  developed  in  this  piogram  for  the  compu- 
tation of  cone  penetration  resistance  on  the  assumption  that  the  soil  above  the  base  of  the 
cone  acts  as  a surcharge  (by  its  weight  only)  as  illustrated  in  Fig.  13.  The  slip  line  field 
shown  in  the  Figure  is  the  solution  of  the  differential  equations  of  plasticity  for  the  axially 
symmetric  case  obtained  by  numerical  solution  methods.  This  reseai'ch  project  revealed 
that  for  certain  cone  angles  and  c - 0 values  it  was  not  possible  to  close  the  slip  line  field 
at  the  apex,  and  for  some  conditions,  overlapping  of  the  slip  lino  field  occurred  (Fig.  14a), 
Such  overlapping  is  physically  not  permissible  since  two  difieront  stress  states  cannot 
exist  in  the  soil  at  the  some  location  and  some  time. 

To  resolve  these  problems  further  theoretical  studies  were  made  for  the  present 
project.  Those  indicated  that  for  the  30®  apex  angle  WES  cone  axial  symmetry  requires 
a different  direction  of  the  major  principal  stress  at  the  apex  than  that  specified  as 
boundary  condition  ior  the  face  of  the  cone,  including  the  apex,  as  the  direction  corre- 
sponding to  the  assumed  interface  friction  angle.  Thus,  in  the  immediate  vicinity  of  the 
apex  plastic  state  and  axial  symmetry  pose  conflicting  requirements  and,  therefore,  some 
other  than  a plastic  state  (probably  rigid)  must  exist  there.  Fortimately,  for  the  total  cone 
resistance  this  problem  has  little  significance  since  the  area  involved  Is  very  small.  In 
the  following  analyses  l/lOO  of  the  base  area  (within  one  tenth  of  the  base  radius)  was 
excluded  from  the  slip  line  field  computations  to  solve  the  problem.  It  is  believed  that 
the  exclusion  of  such  a small  area  does  not  affect  the  results  appreciably. 
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Fig.  13.  Application  of  Plaiticity  Theory  to  Cone  Penetration  Problem.  — 
Weight  of  Soil  Above  Base  Level  of  Cone  Acts  as  Surcharge  (w). 
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1 

'«i 
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•)  OVERLAP  IF  SOIL  ABOVE  BASE  LEVEL  OF  CONE  IS  TREATED  AS  SURCHARGE 


b)  NO  OVERLAP  WITH  INCREMENTAL  PENETRATION  METHOD 


Fifl'  14.  Slip  Un«  Fisldi  for  Con*  Pen«tr<itlon  Problem. 
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The  overl^ping  problem  was  resolved  by  treating  the  depth  effect  differently  than 
in  the  earlier  project.  Instead  of  assuming  that  the  soil  above  the  base  of  the  cone  acts 
by  its  wei^  only,  it  is  assiuned  that  the  stress  state  in  the  soil  caused  by  t>''’  penetration 
of  the  cone  is  "locked  in. " The  penetration  problem  is  solved  in  an  Incremental  way; 
starting  from  the  surface,  the  stress  state  at  a depth  A z is  determined  from  the  slip  line 
field  and  in  the  next  step  the  slip  line  field  is  computed  for  the  condition  that  the  base  of 
the  cone  is  at  depth  A z and  the  boundary  conditions  at  that  level  outside  the  base  area 
(previously  assumed  as  surcharge  corresponding  to  A z depth)  are  determined  from  the 
"locked  in"  stress  state.  The  stress  state  at  this  depth  has  been  determined  in  the  pre- 
vious slip  line  field  computation.  When  the  boundary  conditions  correspond  to  this  "locked 
in"  stress  state,  there  is  no  overtyping  (Fig.  14b).  Variation  of  cone  penetration  resist- 
ance with  depth  computed  by  this  method  is  shown  in  Fig.  15  for  three  values  of  the  co- 
hesion. 

This  method  is  suitable  for  determination  of  the  variation  of  cone  penetration 
resistance  with  depth,  if  the  Coulomb  strength  parameters,  the  unit  weight  ( 7),  and  the 
interface  friction  angle  ( (p  ) are  known.  It  is  Interesting  to  note  that  even  for  a.  p = 20° 
material,  the  cone  penetration  resistance  remains  constant  beneath  a certain  depth.  This 
eylalns  the  often  observed  phenomenon  that  cone  index  signatures  in  frictional-cohesive 
soils  often  resemble  those  obtained  in  purely  cohesive  soils. 

This  method  may  be  used  to  establish  Coulomb  strength  parameters  from  cone  index 
profiles  by  trial  and  error.  Further  theoretical  and  experimental  research  (including 
interface  friction  angle  measurements)  is  needed  to  develop  a procedure  that  uniquely 
converts  cone  index  profiles  to  Coulomb  strength  parameters.  It  is  strongly  recommended 
that  further  field  cone  penetration  tests  be  conducted  in  such  a way  that  the  complete  depth 
profile  of  cone  resistance  be  available  for  possible  conversion  to  Coulomb  strength 
parameters. 

EFFECT  OF  THE  PASSAGE  OF  LEAD  WHEEL  ON  THE  PROPERTIES  OF  SOIL 
ENCOUNTERED  SY  TRAILING  WHEEL 

Frictional  Soils 

The  passage  of  a tire  generally  increases  the  relative  density.  This  increase  in 
relative  density  depends  on  the  limit  pressure.  The  estimated  relative  density  after 
passage  of  the  lead  whesl  is 

D'^.  = P^  (l-D^)/50  (5) 


33 


where  = initial  relative  density 

D'^  = relative  density  after  passage  of  lead  wheel 
Pj  = limit  pressure  for  lead  wheel 

The  maximum  value  of  Dj,  computed  by  Eq.  (5)  is  set  at  = 1. 1 (D^,  = 1 is  the 
maximum  relative  density  obtainable  by  laboratory  procedures.  In  the  field  somewhat 
higher  densities  may  be  obtained).  ^ 

Cohesive  Solis 

Soils  in  this  category  are  close  to  100%  saturation  at  which  preloading  by  a lead 
tire  does  not  increase  the  soil  strength.  Multiple  passage  of  vehicles  may  even  destroy 
some  of  the  structural  strength  of  this  type  of  soil. 

Cohesive- Frictional  (c-^)  Soils 

Cohesive-frictional  soils  encountered  at  the  surface  are  generally  not  saturated. 
The  presence  of  air  voids  in  these  soils  allows  the  soil  to  compact  under  the  stresses 
applied  by  the  lead  wheel  to  the  soil.  If  the  degree  of  saturation  is  less  than,  say,  about 
85%  then  it  may  be  assumed  that  the  stresses  in  the  soil  are  '’effective  stresses",  l.e. , 
no  part  of  the  applied  stresses  is  carried  by  pore  water  or  air  pressure.  Once  the  soil 
acquired  some  strength  in  its  effective  stress  history,  a major  portion  of  this  strength 
remains  "locked  in"  oven  after  the  stresses  are  released.  This  concept  may  bo  used  for 
the  estimation  of  the  increase'  in  soil  strength  duo  to  the  action  of  the  lending  wheel 
(Fig.  16).  The  maximum  normal  stress  that  the  soil  experiences  beneath  the  lead  wheel 
dotermlnea  the  strength  that  remains  partially  locked  in  the  soil.  The  degree  of  "locking 
in"  is  expressed  by  the  coefficient  K that  is  applied  to  the  friction  angle  in  the  stress 

a 

range  up  to  the  estimated  max.  normal  stress.  IC  may  be  determined  by  triaxial  tests 

s 

that  duplicate  the  stress  path  in  the  soil  duriiig  passage  of  the  load  wheel. 

Another  way  of  obtaining  information  on  the  effect  of  compaction  on  the  strength 
properties  of  soil  would  be  to  make  cone  ptaiotration  teats  in  the  rut  of  vehicles.  In  the 
f : dd  there  is  generally  some  off-road  vehicle  that  carries  t'iie  crew.  It  would  require 
little  additional  work,  to  make  cone  ixjnotratijn  tests  in  the  rut  of  the  vehicle  and  record 
it  togotlwr  with  the  Information  on  the  vehicle  tiro  characteristics. 
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Fly.  16.  Eitiimtlon  of  th«  Eff«ot  of  Cotn)Motion  by  th«  Load  Whaul  on  tha 
Strangth  Propartiai  of  Soil  bnoounurad  by  tha  Trailing  Whaal. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 


A vehlole  porformance  model  has  been  developed  that  Incorporates  all  essential 
interactive  features  of  vehlole  soil  interaction  that  affect  vehicle  peirformanoca  The  model 
simulates  vehlole  porformance  with  sufficient  accuracy  lor  both  periormonce  prediction 
and  parametric  analyses  of  conceptual  and  existing  vehicles.  A braked  tlre-soll  inter- 
action model  has  also  been  developed  to  complement  the  vehicle  performance  model. 

It  is  recommended  that  the  Interactive  oapobllity  of  the  vehicle  performance 
model  be  used  to  establish  applied  torque-motion  resistance  relationships  for  the  various 
military  vehicles.  The  use  of  such  relationships  in  the  AMC  '74  Mobility  Model  would  en- 
hance the  performance  simulation  therein  without  significant  chan.'^es  in  the  present  struc- 
ture of  the  model.  It  is  also  recommended  that  the  present  4x4  vehicle  performance  mod- 
el be  expanded  to  multiaxle  configurations. 

In  the  soil  categories  presently  used  In  the  AMC  '74  Mobility  Model  frictional- 
cohesive  soils  are  not  included.  It  is  recommended  that  a new  soil  category  comprising 
these  soils  be  included  in  the  Model  and  the  vehicle  performance  model  be  used  for  the 
prediction  of  tractive  performance  in  this  soil  category.  For  the  characterization  of  soils 
in  this  category  it  is  necessary  to  establish  their  Coulomb  strength  parameters.  For  this 
purpose  it  is  recommended  that 

a)  Complete  penetration  resistance-depth  profiles  be  obtained  in  future  field 
investigations,  both  on  virgin  soil  and  after  the  passage  of  a specified  vehicle 

b)  I'urther  research,  based  on  the  new  method  of  analysis  of  continuous  cone 
penetration  (reported  in  Section  4),  be  conducted  to  develop  methods  for  the  direct  con- 
version of  resistance  profiles  to  Coulomb  strength  parameters. 
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